Q-scan-analysis of the neutron scattering in iron-based superconductors 
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We propose a way to determine the pairing state of the iron pnictide superconductors exploiting 
the momentum (Q) scan of the neutron scattering data. We investigate the spin susceptibility 
in the s± and s++ superconducting states for various doping levels using the effective five-orbital 
model and considering the quasiparticle damping. The peak position of the intensity shifts from 
the position on the line Q x = n to that on the line Q y = as the doping level is decreased from 
electron doping to hole doping. We find that the Q-dependence of the ratio of the intensity in the 
superconducting state to that in the normal state is qualitatively different between the s±-wave and 
s++-wave pairings. We propose to investigate experimentally this ratio in Q-space to distinguish 
the two pairing states. 

PACS numbers: 74.20. Rp, 78.70.Nx, 74.70.Xa 



I. INTRODUCTION 

The discovery of the iron-based superconductors has 
attracted considerable attention because of the high su- 
perconducting transition temperature^. The possibil- 
ity of a peculiar unconventional pairing state has also 
been an issue of great interest. In fact, spin-fluctuation- 
mediated s±-wave pairing state has been proposed at the 
early stage of the study 2 ^— , where the superconducting 
gap is fully open, but changes its sign across the wave 
vector that bridges the disconnected Fermi surfaces. The 
sign change occurs because repulsive pairing interaction 
arises from the spin fluctuations that develop around the 
Fermi surface nesting vector. 

There are many experimental results which suggest 
that the order parameter is fully gapped in a num- 
ber of iron-based superconductors, such as the penetra- 
tion depth measurements 8 - - — , the angle resolved photo- 
emission spectroscopy^—, and the scanning tunneling 
microscopy/spectroscopy (STM/STS) 1 ^. There are also 
experimental suggestions that the superconducting gap 
has an unconventional form with sign change. The nu- 
clear magnetic relaxation rate lacks the coherence peak 
below T,~~—, which suggests that there is a sign change 
in the order parameter. In addition, integer and half- 
integer flux-quantum transitions in composite niobium- 
ion pnictide loops have also been observed 1 ^, which again 
suggests the presence of the sign change in the order pa- 
rameter. In Ref. [Til STM/STS measurements have been 
performed to detect the quasiparticle interference origi- 
nating from the s±-gap. These experiments seem to be 
consistent with the s±-wave pairing scenario. 

On the other hand, there has been a debate concern- 
ing the sensitivity of T c against impurities. There have 
been a number of studies regarding the issue of whether 
the pairing state in the iron pnictides is robust against 
impurities or not^Sr— . It has been pointed out that the 
suppression of T c by impurities is too week for a pair- 



ing state with a sign change in the gap in some papers 
e.g., Ref. [23|. A calculation based on a five-band model 
by Onari et al. has supported this theoretically 2 ^, al- 
though the strength of the impurity potential adopted 
there is large compared to those calculated from first 
principles 2 ^. As a possible pairing state that is robust 
agains impurities, the so-called s ++ -state, where the gap 
does not change its sign between the Fermi surfaces, has 
been proposed 2 ^—. 

As for the probe to determine the pairing state, it has 
been proposed at the early stage d 30 ' 31 that the observa- 
tion of neutron scattering resonance at the nesting vector 
of the electron and hole Fermi surfaces is one of the use- 
ful ways to determine whether there is sign change of the 
gap between these disconnected Fermi surfaces. In fact, 
neutron scattering experiments have observed a peak-like 
structure in the superconducting stated— . This has 
been taken as strong evidence for the sign change in the 
superconducting gap. On the other hand, Onari et al. 
later took into account the quasiparticle damping effect 
in the calculation of the dynamical spin susceptibility, 
and showed that a peak-like enhancement over the nor- 
mal state values can be obtained even in the s++ state, 
which is due to the suppression of the normal state sus- 
ceptibility originating from the damping 38,39 . In these 
studies, they claim that the enhancement that originates 
from the resonance of the s±-wave pairing is too strong 
to explain the experimental observations, and the 'hump- 
like structure' of the intensity in the s ++ -wave pairing is 
more likely to be the origin. However, a quantitative 
comparison between theories and experiments should ac- 
tually be very difficult since there are many parameters 
that can affect the intensity, e.g., the doping- level, the 
band-structure, the strength of the electron interaction, 
and among all, the adopted theoretical method (random 
phase approximation, fluctuation exchange, etc.) There- 
fore, we need to look for a more qualitative difference 
between the peak structures in the two pairing states, 



which can be understood intuitively. Recently, we have 
proposed to investigate experimentally the wave vector 
~ (71", 7r) in the unfolded Brillouin zone, in addition to 
the usually considered ~ (it, 0)/(0, 7r), to distinguish the 
two pairing states qualitatively^. In this paper, we ex- 
tend this study, and obtain a Q (momentum)-.E(energy) 
scan of the dynamical susceptibility. The momentum is 
scanned along the whole symmetric lines. We also in- 
vestigate the doping concentration dependence, where a 
recent experimental as well as theoretical studies in the 
normal state have revealed an electron-hole asymmetry of 
the incommensurability of the spin fluctuatio n 44 ^ 5 . We 
propose that by looking at the superconducting to nor- 
mal state ratio of the dynamical susceptibility, the two 
states give qualitative difference reflecting the difference 
in the mechanism of the occurrence of the peak like struc- 
ture. The ratio is barely momentum dependent in the 
s++ state, while in the s± state it is maximized around 
the wave vector at which the normal spin susceptibility 
is maximized. 

This paper is organized as follows. In Sec. II, four 
effective models for electron- or hole- doped systems and 
the formulation of the spin susceptibility are presented. 
We introduce the multi-orbital random phase approxi- 
mation (RPA) with the quasiparticle damping. In Sec. 
Ill, the calculation results are presented. We systemat- 
ically investigate the doping-level and scattering vector 
dependences, and point out a qualitative difference be- 
tween s±-wave and s++-wave states. In Sec. IV. the 
conclusion is given. 



II. FORMULATION 
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FIG. 1. Fermi surfaces for (a) heavily electron-doped (n ~ 
6.3), (b) optimally electron-doped (n ~ 6.1), (c) lightly 
electron-doped (n ~ 6.04), and (d) optimally hole-doped 
(n ~ 5.8) cases. 



B. Spin susceptibility 

The dynamical spin susceptibility is given by 

x(Q,e) = Y / xTc(Q,e). (1) 



A. The effective models 

We introduce two-dimensional five-orbital models of 
1111 materials obtained in the unfolded Brillouin zone 7 , 
where the x- and y- axes are taken in the Fe-Fe bond 
direction. To investigate the doping-level dependence, 
we consider four situations : lightly electron-doped (the 
band filling n ~ 6.04) , optimally electron-doped (n ~ 
6.1), hcavily-electron doped (the band filling is n ~ 6.3), 
optimally hole-doped (n ~ 5.8) cases. For the electron 
doped cases, we use the model of LaFeAsO, while for 
the hole doped case, to mimic the band structure of the 
122 materials, where the xy hole Fermi surface around 
the wave vector (jr, tt) in the unfolded Brillouin zone is 
robust, we use a two dimensional model of NdFeAsO 43 . 
This is just for simplicity in the calculation, namely, the 
model of the 122 systems have strong three dimension- 
ality, which is difficult to cope with within the present 
formalism. We expect that the three dimensionality it- 
self does not affect the conclusion of the present study. 
The Fermi surfaces of each dopings are shown in Fig. [T] 



Here, XcciQ'E) denotes an orbital-dependent spin sus- 
ceptibility. The inelastic neutron scattering intensity is 
proportional to x"(Qi E) which is the imaginary part of 
x(Q, E). In the multi-orbital RPA, an orbital-dependent 
spin susceptibility XcciQ^E) is written as 

-i aa 

X™{Q,E)= (l-MQ^)U a )- 1 ico{Q,E) , (2) 

J cc 

where xo{QtE) is the bare spin susceptibility expressed 
as 



lxo(Q, E)C = -J2H [K^ik, k + Q)xZ(k, k + Q,E) 

k vv' 

+MX'J (fc, k + Q)x^i (k, k + Q, E) . 

(3) 

Here, Xgg ( . F ^(fe,fc + Q,E) denotes the normal (anoma- 
lous) part of the band-dependent BCS spin susceptibility 
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at zero-temperature (E > 0) with E% — y/e 
K\ 2 = (1 + ^/^fe)/2, K\ 2 = (1 - c£/^)/2, «* = 
A^/ and e£ is the i/-th band energy measured rela- 
tive to the Fermi energy. M^' c f{k,k+Q) {M^'J(k,k + 
Q)) is given by 

M a t d G (fc, k + Q) = u: v {k)U bu> {k + Q)U* c Ak + Q)U dv (k), 

(6) 



MXT (fc, fe + Q) = U* ay {k)U bu ,{k + Q)U clJ {W*dv>{k + Q), 

(7) 

with the unitary matrix U (fe) which diagonalizes the 
Hamiltonian in the orbital basis. Here, we introduce 
the band-index v whose energy satisfies the relation 



[y > v'). For the 's++-wave', we take A 2 



A a = A 4 = A and for 's±-wave' A 2 = A a = -A 4 = A . 
As done in Ref. HH, we introduce a Gaussian cutoff for 
the gap A^ = A u exp{-[e^/A^] 2 }, and take AE = 4A . 

We employ the orbital dependent interactions in the 
form U s = af/Miyake- Here, C/Miyake are the orbital depen- 
dent interactions obtained from the first-principles calcu- 
lation by Miyake et al£^ In the RPA calculation, since 
realistic values of the interaction results in very large spin 
fluctuations, we multiply all the electron-electron interac- 
tion by a reduction factor a. We set a — 0.5 for n ~ 6.1, 
a = 0.45 for n ~ 6.04, a = 0.55 for n - 6.3, and a = 0.4 
for n ~ 5.8. These values of a are close to the maxi- 
mum value that can be adopted in the present formalism 
(larger values give magnetic transition), which is depen- 
dent on the band filling. The variance of the maximum 
value of a is actually largely due to the approximation 
adopted here (RPA), where the self energy corrections 
are not taken into account. Nevertheless, we have con- 
firmed that our results do not change qualitatively if we 
take other values of the reduction factor. 

We consider the quasiparticle damping rj^ in the form 



v k.a = m ax(7fc,7fc +Q ), 



(8) 
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This form of the damping was introduced in Ref. |38j from 
the requirement that the damping should not be present 
in the low energy regime of less than 3A, which is equal 
to the particle-hole excitation gap plus one-particle gap. 



In the high energy regime, the damping should be es- 
sentially the same as that in the normal state, and the 
present form interpolates the low and high energy regimes 
in a simple manner. To be more strict, the damping of 
the superconducting state in the high energy regime can 
be different from that in the normal state above T c be- 
cause the temperature is different. However, it is known 
experimentally that the dynamical susceptibility in the 
normal and the superconducting states nearly coincide 
with each other in the high energy regime^, indicating 
that the temperature dependence of the damping in the 
high energy regime is not so strong. Therefore, we take 
the same value of 70 for the normal and the supercon- 
ducting states. Although we use this specific form of the 
damping in the actual calculation, the present results are 
not qualitatively dependent on the details of the damp- 
ing. 



III. RESULTS 

We investigate the Q-E map of the spin susceptibility 
x"(Q, E). Since the commensurate vector is Q — (n, 0), 
we calculate Q-dependence along two lines in Q-map. 
One is the line on the Q^-axis, i.e., Q = (Q x , 0), and the 
other is the line with Q x — n, Q — (tt, Q y )- In this paper, 
the smearing factor and the magnitude of the supercon- 
ducting gap are taken to be 77 = 0.5 meV, A = 10 meV, 
respectively. We take the quasiparticle damping in the 
normal state as 70 = 10 meV, which was estimated from 
the temperature dependence of the resistivity observed 
experimentally^. To cope with the realistic magnitude 
of the superconducting gap, we take 4096 x 4096 fe-point 
meshes throughout the paper. 

First, we will show the results in the normal state to see 
the peak position of the spin susceptibility. The incom- 
mensurability is found to be asymmetric between electron 
and hole doped cases, as was found experimentally^ 4 - as 
well as theoreticall y 45 ' 46 . Then, we will present the re- 
sults in the superconducting state to investigate the dif- 
ference between s-t-wave and s++-wave states in Q-maps. 
We find that the Q-map of the superconducting/normal 
ratio is qualitatively different between the two states. 



A. Normal state 

We show the results in the normal state in Fig. O We 
find that the peak position of the intensity shifts from the 
position on the line Q x = tt to that on the line Q y = as 
we decrease the electron doping level, and go into the hole 
doping regime (see Fig. [3]). This electron- hole asymme- 
try of the incommensurability of the spin fluctuations has 
been attributed to the multiorbital nature of the Fermi 
surfac o 44 ' 45 . It should also be noted in the heavily elec- 
tron doped regime of n — 6.3, the incommensurability 
is very large, coming somewhat close to (tt/2,tt). This 
structure comes partially from the interaction between 
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FIG. 2. Q-i? map of the spin susceptibility x"(Q> E) m the 
normal state for heavily electron-doped (the band filling is 
n ~ 6.3), optimally electron-doped (n ~ 6.1), lightly electron- 
doped (n ~ 6.04), and optimally hole-doped (n ~ 5.8) cases. 
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FIG. 4. Q-i? map of the spin susceptibility x"(Q> E) m the 
s±-wave pairing state for the four band filling cases. 
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FIG. 3. Schematic figure of the peak-position shift in Q-space 
with decreasing the doping level. 



two electron pocketsZ^ 3 -, but there turns out to be also 
a contribution from the interaction between the electron 
and the barely present hole Fermi surfaces (as we shall 
see more clearly in the superconducting state) , which are 
not well nested at this doping level. In this sense, we are 
using the term "nesting" in a rather broad sense of the 
term. Another point to be noted in the heavily electron 
doped regime (n = 6.3) is the absence of intensity at low 
energies, which reflects the fact that the hole Fermi sur- 
faces are barely present. As for the interaction strength 
dependence, we have confirmed that the peak position 
does not depend on the interaction reducing factor a. 



electron-doped case (n ~ 6.1) gives the largest peak in- 
tensity among these four doping levels. Therefore, we 
point out once again that it is difficult to determine the 
pairing state from quantitative comparison between ex- 
periments and theory since the peak intensity depends 
on the doping level n as well as the strength of the elec- 
tron interaction U s . Another point that should be noted 
here is that the resonance peak exists in the heavily elec- 
tron doped regime n = 6.3, reflecting the sign change of 
the gap between electron and hole Fermi surfaces. This 
means that there is a certain contribution to the structure 
close to (7r,7r/2) also from the "badly nested" electron 
and hole Fermi surfaces, as mentioned in the section for 
the normal state. A resonance at a wave vector close to 
("7r, 7r/2) in heavily overdoped systems, as is observed for 
instance in Ref.— , is interpreted as a signature of d-wave 
pairing 4 ™ where the gap sign changes between electron 
pockets, but the present study suggests that it may also 
be interpreted in terms of s± pairing originating from 
badly nested electron and hole pockets^. 



B. Superconducting state: s±-wave pairing 



C. Superconducting state: s++-wave pairing 



Let us now consider the superconducting state with 
s±-wave pairing. We show the Q-E maps of the spin 
susceptibility %"(Q, i?) in Fig. HI We can clearly see 
the shift of the peak position with decreasing the doping 
level, i.e., the resonance peak moves in accord with the 
normal state peak position. It should be noted that the 
peak intensity also depends on the doping level as well 
as the interaction strength. For example, the peak at the 
incommensurate vector Q — (it, 7r/8) for the optimally 



Next, we consider the s++-wave state. We show the 
Q-E maps of the spin susceptibility x"{Q,E) m Figs. [5] 
As pointed out previousl y 38 ! , the quasiparticle damping 
makes the hump structure in the s++-wave pairing state. 
The hump structure moves as the doping level is varied, 
which at first sight may look similar to that in the s± 
state. However, there is a large difference between these 
two states, which is revealed by taking the ratio between 
the normal and the superconducting states. 
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FIG. 5. Q-E map of the spin susceptibility x'iQ^E) in the 
s++-wave pairing state for the four band filling cases. 



D. Superconducting to normal ratio 

We now show the Q-E map of the ratio of the spin 
susceptibility in the superconducting state to that in the 
normal state in Figs. [51 As shown in the left panels in 
these figures, the ratio becomes large in a localized regime 
around Q ~ (ir, 0) in the case of the s±-wave pairing. To 
be more precise, one can find in Fig. [7] that the peak 
position of the ratio in the low energy region (E < 2Aq) 
shifts with decreasing the doping level, in accord with 
the peak position of the spin susceptibility in the normal 
state as was shown in Fig. [3] By contrast, in the case 
of ,s ++ -wave pairing, we find that the ratio is weakly 
dependent on Q as shown in the lower panels in Fig. [6] 
This is a qualitative difference between the s±-wave and 
,s ++ -wave pairing. 

The above difference originates from the difference of 
the mechanism of the enhancement between the two 
cases. In the s±-wave pairing state, the enhancement 
originates from the sign change of the gap across the wave 
vector that bridges the Fermi surface, so that the maxi- 
mum value of the ratio in Q-space is localized around this 
wave vector. On the other hand, in the s+_|--wave pair- 
ing state, the enhancement comes from the quasiparticle 
damping effect, which is not momentum dependent. As 
pointed out in previous studie s 39 ' 40 , the dissipation in- 
creases in the high energy region (E > 3Ao) is the origin 
of the decrease of the intensity in the high energy re- 
gion, so that one can see the hump structure as shown in 
Fig.[5]in the mid-range energy region (2Ao < E < 3Ao). 
This means that the enhancement (the superconducting 
to normal ratio) in the s++-wave state and the nesting 
vector are not related. This was the original motivation 
for the proposal to investigate experimentally the wave 
vector ~ (7r,7r) to distinguish between two cases^ -. In 
the present study, it can be more clearly seen that the 
ratio of the enhancement in the s++-wave state depends 
weakly on the wave vector. In addition, we also point out 
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FIG. 6. Q-E map of the ratio of the intensity in the super- 
conducting state to that in the normal state for the four band 
filling cases. Left panel: the result of the s±-wave pairing 
state. Right panel: the result of the s++-wave pairing state. 



that this ratio does not depend on the doping level. This 
result is easily understood by the nature of the quasi- 
particle damping. Thus, if the experimentally observed 
peak structure is coming from the quasiparticle damping, 
the enhancement ratio of the intensity in the supercon- 
ducting state to that in normal state has to be weakly 
momentum dependent. This conclusion is valid in vari- 
ous kinds of the iron-based superconductors. 



IV. CONCLUSION 

In conclusion, we have obtained the Q-map of the spin 
susceptibility in the s± and 3+4. superconducting states 
of the iron-based superconductors by applying multior- 
bital RPA to the effective five-orbital models and con- 
sidering the quasiparticle damping. We have shown that 
the peak position of the intensity shifts from the posi- 
tion on the line Q x = n to that on the line Q y = 
with decreasing the doping level. Considering the diffi- 
culty in quantitatively comparing the theories with the 
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FIG. 7. Q-dependence of the x"{Q> E) in the s±-wave pairing 
state. Left panel: Qz-dependence Q = (Q x ,0). Right panel: 
Q M -dependence Q — (ir,Q v ). 



experiments, we have proposed that the comparison of 
the Q-dependence of the superconducting to normal ra- 
tio of the spin susceptibility is useful in determining the 
pairing state. In the s±-wave state, the ratio is peaked 
around the nesting vector, and it is dependent on the 
doping level, while in the s++-wave pairing state, the ra- 
tio is weakly dependent on the wave vector. Exploiting 
the present proposal, we expect that one can distinguish 
the two pairing states. 



ACKNOWLEDGMENT 



We thank M. Machida, N. Nakai, Y. Ota, for 
helpful discussions and comments. The calculations 
have been performed using the supercomputing sys- 
tem PRIMERGY BX900 at the Japan Atomic Energy 
Agency. 



1 Y. Kamihara, T. Watanabe, M. Hirano, and H. Hosono, J. 
Am. Chem. Soc. bf 130, 3296 (2008). 

2 Y. Bang, and H.-Y. Choi, Phys. Rev. B, 78, 134523 (2008). 

3 K. Seo, B. A. Bernevig, and J. Hu, Phys. Rev. Lett. 101, 
206404 (2008). 

4 M. M. Parish, J. Hu, and B. A. Bernevig, Phys. Rev. B 
78, 144514 (2008). 

5 M. M. Korshunov and I. Eremin, Phys. Rev. B 78, 
140509(R) (2008). 

6 I. I. Mazin, D. J. Singh, M. D. Johannes, and M. H. Du, 
Phys. Rev. Lett. 101, 057003 (2008). 

7 K. Kuroki, R. Arita, H. Usui, Y. Tanaka, H. Kontani, and 
H. Aoki, Phys. Rev. Lett. 101, 087004 (2008). 

8 K. Hashimoto, T. Shibauchi, T. Kato, K. Ikada, R. 
Okazaki, H. Shishido, M. Ishikado, H. Kito, A. Iyo, H. 
Eisaki, S. Shamoto, and Y. Matsuda, Phys. Rev. Lett. 102, 
017002 (2009). 

9 H. Luetkens, H.-H. Klauss, R. Khasanov, A. Amato, R. 
Klingeler, I. Hellmann, N. Leps, A. Kondrat, C. Hess, A. 
Kohler, G. Behr, J. Werner, and B. Biichner, Phys. Rev. 
Lett. 101, 097009 (2008). 

10 L. Malone, J. D. Fletcher, A. Serafin, A. Carrington, N.D. 
Zhigadlo, Z. Bukowski, S. Katrych, and J. Karpinski, Phys. 
Rev. B 79, 140501 (R) (2009). 

11 H. Ding, P. Richard, K. Nakayama, T. Sugawara, T. 
Arakane, Y. Sekiba, A. Takayama, S. Souma, T. Sato, T. 
Takahashi, Z. Wang, X. Dai, Z. Fang, G. F. Chen, J. L. 



Luo, and N. L. Wang, EuroPhys. Lett. 83, 47001 (2008). 

12 K. Nakamura, T. Sato, P. Richard, Y.-M. Xu, Y. Sekiba, 
S. Souma, G. F. Chen, J. L. Luo, N. L. Wang, H. Ding, 
and T. Takahashi, Europhys. Lett. 85 67002 (2009). 

13 D. V. Evtushinsky, D. S. Inosov, V. B. Zabolotnyy, 
A.Koitzsch, M. Knupfer, B. Buchner, M. S. Viazovska, G. 
L. Sun, V. Hinkov, A. V. Boris, C. T. Lin, B. Keimer, A. 
Varykhalov, A. A. Kordyuk, and S. V. Borisenko, Phys. 
Rev. B 79, 054517 (2009). 

14 T. Hanaguri, S. Nittaka, K. Kuroki, and H. Takagi, Science 
328, 474 (2010). 

Y. Nakai, K. Ishida, Y. Kamihara, M. Hirano, and H. 
Hosono, J. Phys. Soc. Japan, 77 073701 (2008). 

16 H.-J. Grafe, D. Paar, G. Lang, N. J. Curro, G. Behr, J. 
Werner, J. Hamann-Borrero, C. Hess, N. Leps, R. Klin- 
geler, and B. Buchner, Phys. Rev. Lett. 101, 047003 
(2008). 

17 H. Mukuda, N. Terasaki, H. Kinouchi, M. Yashima, Y. 
Kitaoka, S. Suzuki, S. Miyasaka, S. Tajima, K. Miyazawa, 
P. Shirage, H. Kito, H. Eisaki, and A. Iyo, J. Phys. Soc. 
Jpn. 77, 093704 (2008). 

18 N. Terasaki, H. Mukuda, M. Yashima, Y. Kitaoka, K. 
Miyazawa, P. M. Shirage, H. Kito, H. Eisaki, and A. Iyo, 
J. Phys. Soc. Jpn. 78, 013701 (2009). 

19 C.-T. Chen, C. C. Tsuei, M. B. Ketchen, Z.-A. Ren, and 
Z. X. Zhao, Nature Physics 6, 260 (2010). 



7 



20 Y. F. Guo, Y. G. Shi, S. Yu, A. A. Belik, Y. Matsushita, M. 
Tanaka, Y. Katsuya, K. Kobayashi, I. Nowik, I. Felner, V. 
P. S. Awana, K. Yamaura, and E. Takayama-Muromachi, 
Phys. Rev. B 82, 054506 (2010). 

21 Y. Nakajima, T. Taen, Y. Tsuchiya, T. Tamegai, H. Ki- 
tamura, and T. Murakami, Phys. Rev. B 82, 220504(R) 
(2010). 

22 S. Kitagawa, Y. Nakai, T. Iye, K. Ishida, Y. F. Guo, Y. 
G. Shi, K. Yamaura, and E. Takayama-Muromachi, Phys. 
Rev. B 83, 180501(R) (2011). 

23 M. Sato, Y. Kobayashi, S. C. Lee, H. Takahashi, E. Satomi, 
and Y. Miura, J. Phys. Soc. Jpn. 79 (2010) 014710. 

24 Y. K. Li, X. Lin, Q. Tao, C. Wang, T. Zhou, L. J. Li, Q. 
B. Wang, M. He, G. H. Cao, and Z. A. Xu, New J. Phys. 
11 053008 (2009). 

25 S. Onari and H. Kontani, Phys. Rev. Lett. 103 (2009) 
177001. 

26 K. Nakamura, R. Arita, and H. Ikeda, Phys. Rev. B 83, 
144512 (2011). 

27 H. Kontani and S. Onari, Phys. Rev. Lett. 104, 157001 
(2010). 

28 Y. Yanagi, Y. Yamakawa, N. Adachi, and Y. Ono, J. Phys. 
Soc. Jpn. 79,123707 (2010). 

29 T. Saito, S. Onari, and H. Kontani, Phys. Rev. B 82 
144510 (2010). 

30 TA. Maier and D. J. Scalapino, Phys. Rev. B 78, 
020514(R) (2008). 

31 T. A. Maier, S. Graser, D. J. Scalapino, and P. Hirschfeld, 
Phys. Rev. B 79, 134520 (2009). 

32 A. D. Christianson, E. A. Goremychkin, R. Osborn, S. 
Rosenkranz, M. D. Lumsden, C. D. Malliakas, I. S. 
Todorov, H. Claus, D. Y. Chung, M. G. Kanatzidis, R. 
I. Bewley and T. Guidi, Nature 456 930 (2008). 

33 Y. Qiu, W. Bao, Y. Zhao, C. Broholm, V. Stanev, Z. 
Tesanovic, Y. C. Gasparovic, S. Chang, J. Hu, B. Qian, M. 
Fang, and Z. Mao, Phys. Rev. Lett. 103 067008 (2009). 

34 D. S. Inosov, J. T. Park, P. Bourges, D. L. Sun, Y. Sidis, A. 
Schneidewind, K. Hradil, D. Haug, C. T. Lin, B. Keimer 
and V. Hinkov, Nature Phys. 6 178 (2010). 

35 J. Zhao, L.-P. Regnault, C. Zhang, M. Wang, Z. Li, F. 
Zhou, Z. Zhao, C. Fang, J. Hu and P. Dai, Phys. Rev. B 



81, 180505(R) (2010). 

36 M. Ishikado, Y. Nagai, K. Kodama, R. Kajimoto, M. 
Nakamura, Y. Inamura, S. Wakimoto, H. Nakamura, M. 
Machida, K. Suzuki, H. Usui, K. Kuroki, A. Iyo, H. Eisaki, 
M. Arai, S. Shamoto, Phys. Rev. B 84, (2011) 144517. 

37 C. Zhang, M. Wang, H. Luo, M. Wang, M. Liu, J. Zhao, 
D. L. Abernathy, T. A. Maier, K. Marty, M. D. Lumsden, 
S. Chi, S. Chang, J. A. Rodriguez-Rivera, J. W. Lynn, 
Tao Xiang, Jiangping Hu, P. Dai, Scientific Reports 1, 115 
(2011). 

38 S. Onari, H. Kontani, and M. Sato, Phys. Rev. B 81 
060504(R) (2010). 

39 S. Onari, H. Kontani, Phys. Rev. B 84, 144518 (2011). 

40 Y. Nagai, and K. Kuroki, Phys. Rev. B 83, 220516(R) 
(2011). 

41 Y. Nagai, and K. Kuroki. larXiv:1106 .2376 (unpublished); 
S. Onari, and H. Kontani, arXiv: 1107.0748 (unpublished). 

42 T. Miyake, K. Nakamura, R. Arita, and M. Imada, J. Phys. 
Soc. Jpn. 79 (2010), 044705. 

43 K. Kuroki, H. Usui, S. Onari, R. Arita, and H. Aoki, Phys. 
Rev. B 79, 224511 (2009). 

44 C. H. Lee, K. Kihou, H. Kawano-Furukawa, T. Saito, A. 
Iyo, H. Eisaki, H. Fukazawa, Y. Kohori, K. Suzuki, H. 
Usui, K. Kuroki, and K. Yamada, Phys. Rev. Lett. 106, 
067003 (2011). 

45 K. Suzuki, H. Usui, and K. Kuroki, Phys. Rev. B 84, 
144514 (2011). 

46 J.T. Park et al, Phys. Rev. B 82, 134503 (2010). 

47 J.T. Park et al, Phys. Rev. Lett. 107, 177005 (2011). 

48 T. Maier et al, Phys. Rev. B 83, 100515(R) (2011). 

49 Bad nesting of the Fermi surface does not necessarily im- 
ply low T c . In fact, the nesting between electron and hole 
Fermi surfaces is not so good in the system studied by K. 
Kuroki and R. Arita in Phys. Rev. B 64, 024501 (2001) (see 
Fig. 4 of this paper), where a possibility of high T c super- 
conductivity originating from disconnected Fermi surfaces 
was proposed. 



